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摘 要 : 总 初级 生产 力 (GPP) 决 定 了 进入 陆地 生态 系统 的 初始 物质 与 能 量 ,但 由 MODIS CPP 产品 获 


BUN CPP 在 地 表 和 覆盖 复杂 的 黑河 流域 却 不 足以 准确 的 反映 生态 系统 物质 与 能 量 的 分 布 。 因 此 ,利用 
MODIS 影像 数据 ASTER GDEM 数据 .30 mm 分辨 率 土地 利用 履 被 数据 和 中 国 区 域 地 面 气象 要 素 驱 动 
数据 ,驱动 VPM 模型 模拟 黑河 流域 2015 年 5~10 月 的 总 初级 生产 力 ,并 据 此 揭示 黑河 流域 CPP 在 生 
长 季 的 时 空格 局 ,时 间 分 辨 率 为 8 d, 空 间 分 辩 率 为 500 m。 研 究 结 果 表明 :VPM 模型 估算 结果 的 精度 
高 于 MODIS GPP 产品 ,判定 系数 增长 了 45.5% ,总 均 方 根 误差 降低 了 57.0% ;黑河 流域 生长 季 累 积 


CPP 总 体 呈 现 出 中 游 最 高 上 游 其 次 .下 游 最 低 的 显著 空间 分 布 梯度 格局 ;全 境 与 局 部 有 植被 覆盖 区 
域 的 日 GPP 均 呈 先 增加 后 减少 倒 U 型 变化 规律 , 且 前 者 在 7 月 下 旬 达 到 最 高 值 ;植被 覆盖 率 极 低 的 


地 表 区 域 生 长 季 内 的 日 CPP 在 基本 稳定 中 上 下 波动 ,稳定 值 处 于 1 gem? - d Hirt, 
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植被 总 初级 生产 力 ( Gross Primary Production , 
GPP) 是 指 在 单位 时 间 和 单位 面积 上 ,绿色 植被 通过 
光合 作用 所 产生 的 有 机 碳 总 量 "" ,在 评估 生态 系统 
平衡 和 植被 固 碳 能 力 的 过 程 中 起 重要 作用 。 传 统 
GPP 的 测量 方法 主要 有 箱 式 法 、 基 于 涡 度 协 方差 技 
术 的 测定 法 以 及 稳定 同位 素 测 定 方法 等 2 EJL 
十 年 来 ,遥感 技术 得 到 了 巨大 的 发 展 ,其 多 时 空 尺 
度 、 低 成 本 ,高 效率 ,高 精度 的 特点 ,使 区 域 和 全 球 尺 
度 上 估算 CPP 成 为 可 能 “i。 目 前 ,国内 外 学 者 开发 
了 许多 基于 遥感 技术 的 CPP 估算 模型 ,如 GLO- 
PEM {RAI VPM 模型 中 和 BEPS 模型 等。 其 
中 , 光 能 利用 率 模型 将 GPP 调控 因子 以 相对 简单 的 
方法 组 合 ,可 直接 利用 遥感 数据 "" ,已 成 为 CPP 估 
算 的 一 个 重要 发 展 方向 。 

黑河 流域 地 处 一 带 一 路 沿线 ,位 于 干旱 半 干 旱 


VPM 模型 ; SRM; 总 初级 生产 力 ; 时 空格 局 ; 黑河 流域 


高 寒 草 旬 和 中 游 玉 米 的 GPP 动态 进行 了 估算 ;LI 
等 9 采用 REG-PEM 模型 对 黑河 流域 临 泽 .盈科 和 
大 野 口 的 CPP 进行 了 估算 ; 闫 敏 等 "中 利用 改进 的 
MOD17 模型 对 黑河 流域 上 游 的 植被 GPP 进行 了 估 
算 。 然 而 ,目前 这 些 研究 多 集中 于 黑河 流域 中 上 游 ， 
且 研 究 区 多 为 涡 动 站 通 量 贡献 区 范围 ,对 于 黑河 流 
域 下 游 和 流域 整体 的 CPP 估算 研究 则 相对 较 少 。 
因此 ,本 文 尝试 对 黑河 流域 全 境 的 CPP 进行 逐 像 元 
的 估算 ,并 据 此 探索 黑河 流域 物质 与 能 量 的 空间 分 
布 与 时 间 动 态 规律 。 

MODI7 产品 基于 PSN 模型 与 MODIS 数据 对 全 
球 的 CPP 进行 了 逐 像 元 的 估算 ,目前 ,国家 、 省 和 地 
区 等 大 区 域 尺度 的 CPP 研究 多 使 用 MOD17 产品 ， 
但 该 产品 的 CPP 估算 精度 并 不 高 ,因此 ,本 文 的 研 
究 旨 在 提供 一 种 适用 于 复杂 地 表 景 观 与 气候 情况 的 


地 区 ,是 典型 的 内 陆 河 流域 ,具有 3 个 不 同 的 环境 单 
元 ,拥有 极 高 的 研究 价值 。 近 年 来 ,国内 外 许多 学 者 
针对 黑河 流域 的 CPP 估算 展开 了 大 量 的 研究 : 
WANG 等 中 利用 EC 和 VPM 模型 对 黑河 流域 上 游 


© ” 收 稿 日 期 : 2019-03-25; 修订 日 期 : 2019 -08 -14 


高 精度 大 区 域 尺度 CPP 估算 方法 。 本 文 以 黑河 流 
域 为 研究 区 ,利用 可 公开 访问 与 获取 的 多 源 数据 驱 
动 VPM 模型 , 估算 2015 年 植被 生长 季 (5 ~ 10 
FA") #9 8 d 平 均 日 CPP, 并 据 此 揭示 其 时 空格 局 ， 
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研究 成 果 对 于 保护 脆弱 的 黑河 流域 生态 环境 ,指导 
黑河 流域 沿线 的 生态 修复 工程 均 具 有 重要 意义 。 


1 研究 区 概况 


黑河 流域 是 我 国 第 二 大 内 陆 流 域 ,总 面积 约 为 
1.289 x10? km? ,平均 海拔 超过 1 200 m, 主 流 长 度 
约 为 821 km ,地 处 西部 地 区 ,位 于 河西 走廊 中 
部 ,其 地 理 位 置 位 于 37°41’ ~42°42'N 和 96"42' ~ 
102°00' 福 之 间 ( 图 1)。 黑 河流 域 地 形 从 南 到 北 显 著 
不 同 , 按 照 不 同 的 自然 ,生态 和 气候 特征 可 划分 为 上 
游 中游 与 下 游 3 个 区 域 。 上 游 为 祁连山 脉 至 营 落 
峡 峡 谷 ,地 貌 为 高 寒 草 甸 覆 盖 的 山地 地 貌 ;中 游 为 营 
落 峡 峡谷 至 正义 峡 峡 谷 ,地 貌 为 灌溉 农业 的 绿洲 地 
貌 ;下 游 为 正义 峡 峡 谷 至 居 延 海 , 地 貌 为 戈壁 沙漠 覆 
盖 荡 漠 地 貌 '" 。 上 游 .中游 和 下 游 地 区 的 年 平均 气 
温 分 别 为 -3.1~3.6 1€ 7.0 -8.2 'C 18.0 ~ 10.0 
C ,年 平均 降水 量 分 别 为 350 ~ 450 mm, 80 ~ 120 
mm 和 40 ~60 mm!!! 。 
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图 1 黑河 流域 的 地 理 位 置 及 地 形 与 通 量 站 的 空间 分 布 


Fig.1 Location of Heihe River Basin and the spatial 


distribution of terrain and flux station 
2 多 源 数 据 与 预 处 理 


2.1 遥感 数据 
中 分 状 率 成 像 光谱 仪 ( MODIS) 是 NASA 发 射 的 
Terra 和 Aqua 两 颗 卫 星 上 的 中 分 辨 率 高 光谱 传 感 
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3E 。 本 文 使 用 的 遥感 数据 为 MOD09A1、 
MOD15A2H 和 MODI7A2H 三 种 产品 ,数据 获取 网 
址 为 : http://www. ladsweb. nascom. nasa. gov/。 其 
中 ,MOD09A1 为 MODIS 前 7 个 波段 反射 率 的 8 d 合 
成 产品 ;MOD15A2H 为 8 d 的 地 表 叶 面积 指数 与 光 
合 有 效 辐射 吸收 比率 产品 ; MODI7A2H 是 8 d 的 总 
初级 生产 力 产品 ,空间 分 辩 率 为 500 m, 时 间 分 辩 率 
为 8 d, 涵 盖 时 间 为 2015 年 5 ~ 10 月 。 应 用 蓝 波段 
(459 ~479 nm) ` 红 波段 (620 ~ 670 nm) 、 近 红外 波 
段 (841 ~ 875 nm) 和 短波 红外 波段 (1 628 ~ 1 652 
nm) 的 反射 率 进 行 增强 植被 指数 EVI 和 陆地 水 分 指 
Be LSWI 的 计算 ,如 公式 (1) ~ (2) Bro 。 


2.5 x (py = D.) 
[p + (6 x pra 77.5 x Pie) +1 | 


EVI = 


(1) 


《pm Dae) 
HRS (Pair + Dui) (2) 

式 中 :pw Pred Dus 和 pu 分 别 是 蓝 波段 、 红 波段 . 近 
红外 波段 和 短波 红外 波段 的 反射 率 。 本 文 利用 EN- 
VI 软件 预 处 理 生成 黑河 流域 生长 季 的 8 d 合成 
EVI LSWI LAI 和 GPP 数据 产品 。 
2.2 DEM 数据 

本 文 使 用 的 DEM 数据 为 黑河 流域 ASTER 
GDEM 数据 集 , 空 间 分 辨 率 为 30 m。 数 据 获取 网 址 
为 :http://westdc. westgis. ac. cn/。ASTER GDEM 是 
NASA 和 METI 联合 发 布 的 全 球 数 字 高 程 数据 产品 ， 
覆盖 全 球 99% 陆地 表面 ,高 程 精度 为 7 ~14 m", 
本 文 基于 ASTER GDEM 数据 使 用 聚合 法 生成 
480 m fil 0. 1° ADEM, 
2.3 土地 利用 覆 被 数据 

本 文 使 用 的 土地 利用 覆 被 数据 来 自 黑河 流域 土 
地 利用 覆 被 数据 集 必 ,空间 分 辨 率 为 30 m, 数 据 
获取 网 址 为 :http ://westde. westgis. ac. cn/。 黑 河流 
域 土地 利用 覆 被 数据 集 具有 较 高 的 总 体 精度 , 且 细 
化 了 农作物 信息 ,还 更 新 了 冰 、 雪 等 地 类 信息 ,时 间 
分 辩 率 为 1 月 ,时 间 跨 度 为 2011 一 2015 年 。 本 文选 
用 2015 年 7 月 的 地 表 类 型 覆盖 数据 反映 黑河 流域 
生长 季 的 土地 利用 履 被 情况 ,并 根据 植被 功能 型 对 
其 进一步 重 分 类 。 考 虑 到 黑河 流域 的 C4 作物 主要 
有 玉米 ,而 C4 草 类 植物 覆盖 度 很 低 ,遥感 分 类 结 
难以 识别 ,因此 将 农作物 分 为 C3 作物 与 C4 作物 ， 
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Fig.2 Land use and cover with a spatial resolution 


of 30 m in Heihe River Basin on July 2015 


将 草地 统一 划分 为 C3 草地 ,如 图 2 所 示 。 
2.4 气象 数据 

本 文 使 用 的 气象 数据 来 自 中 国 区 域 地 面 气象 素 
驱动 数据 集 王 ; ,数据 下 载 地 址 为 :http ://westdc. west- 
gis. ac. cn/。 中 国 区 域 地 面 气象 要 素数 据 集 是 一 套 近 
地 面 气象 与 环境 要 素 再 分 析 数 据 集 , 由 中 国 科学 院 青 
藏 高 原 研究 所 开发 ,空间 分 辨 率 为 0.1" ,时 间 分 辨 率 
分 为 3 上 与 1 d, 包 含 气 温 气压、 空气 比 湿 、 全 风速 、 
地 面向 下 短波 辐射 .地 面向 下 长 波 辐射 和 地 面 降水 率 
7 个 要 素 。 本 文选 择 3 分辨 率 瞬 时 近 地 面 气温 \ 日 
平均 近 地 面 气温 和 日 平均 地 面向 下 短波 辐射 3 个 数 
据 集 ,时 间 跨 度 为 2015 年 5~10 月 。 

使 用 日 最 高 和 日 平均 近 地 面 气温 的 平均 值 作为 
日 间 平 均 近 地 面 气温 ,其 中 取 每 日 3 分辩 率 瞬时 
近 地 面 气温 中 的 最 大 值 作为 日 最 高 近 地 面 气温 。 对 
于 日 平均 地 面向 下 短波 辐射 ,使 用 双 线 性 插值 算法 
将 空间 分 辨 率 从 0.1° 转 换 为 500 m。 对 于 日 间 平 均 
近 地 面 气温 , 除 使 用 双 线 性 插值 算法 外 ,还 需要 考虑 
海拔 的 影响 ,假设 温度 的 下 降 率 为 6. 50 ,进行 高 
程 订 正 。 
2.5 ”地 面 观 测 数 据 

“黑河 实验 ”在 2010 一 2011 年 完成 总 体 设 计 与 
准备 , 自 2012 年 起 开始 进入 持续 观测 时 期 ,2012 年 
依次 在 中 游 、 上 游 和 下 游 进行 了 加 强 试验 ,全 流域 的 


持续 观测 从 2013 年 开始 至 2015 年 结束 '” 。 本 文 
的 地 面 观测 数据 来 自 " 黑 河 生态 水 文 遥 感 试 验 (Hi- 
WATER)”, 由 8 个 观测 站 的 涡 动 相关 仪 测量 ,测量 
时 间 为 2015 年 全 年 ,站 点 空间 分 布 见 图 1。 涡 动 相 
关 仪 的 原始 观测 数据 为 10 Hz, 发 布 的 数据 集 时 间 
步 长 为 30 min, 下 载 地 址 为 http://westd-c. westgis. 
ac. cn。 该 数据 集 已 经 过 野 点 值 剔除 延迟 时 间 校 
IE 二 次 坐标 旋转 .频率 响应 修正 .超声 虚 温 修正 和 
密度 (WPL) 修 正 等 步骤 的 处 理 ”1 。 对 通 量 数据 
进行 处 理 与 计算 , 即 可 获得 站 点 通 量 贡献 区 的 总 初 
级 生产 力 CPP 。 


3 VPM 模型 计算 与 评估 方法 


3.1 模型 简介 

VPM 模型 是 以 遥感 数据 为 驱动 变量 的 一 种 光 
能 利用 率 模型 ,该 模型 基于 MODIS 遥感 数据 和 涡 度 
相关 碳 通 量 观测 数据 演变 产生 ,考虑 植被 叶绿素 所 
吸收 的 光 能 有 效 辐射 , 且 FPAR 通过 EVI 估算 ,有 效 
的 提高 了 GPP 估算 的 精确 度 。 

TE VPM 模型 中 ,模型 利用 光 能 利用 率 (LUE， 
E, ) 和 叶绿素 吸收 的 光合 有 效 辐射 (4P4R。 ) 的 乘积 
来 估算 总 初级 生产 力 (CPP) , WARG) HR : 

GPP =e, x APAR 4, (3) 

模型 利用 叶绿素 吸收 的 光合 有 效 辐射 比例 
( FPAR a )2 fh TE APAR o FPARA FA EVI 的 线性 也 
数 估 计 , 其 中 系数 a 设 为 1.0, 如 公式 (4) ~ (5) 所 


=~. [22] 
N 


不 


o 


APAR „ = FPAR,,, x PAR (4) 
FPAR =a x EVI (5) 


模型 利用 最 大 光 能 利用 率 sw 温度 限制 因子 
Tu 和 水 分 胁迫 因子 下 的 乘积 来 佑 算 光 能 利用 
率 ,如 公式 (6) 所 示 : 


Eg = Emax x Td x W alar (6) 


在 VPM 模型 中 ,sw 是 特定 生物 群落 所 特有 的 
参数 ,一 般 利 用 土地 利用 覆 被 数据 通过 查找 表 获 取 。 
参考 前 人 研究 '””| ,黑河 流域 植被 最 大 光 能 利用 率 
£y WK 1 所 示 。 

在 VPM 模型 中 ,温度 限制 因子 7 由 陆地 生 
态 系统 模型 (TEM ) 中 使 用 的 公式 (7) 887 : 


( T B Tu ) x ( T 7 rm ) 
T. alar 一 (7 ) 
[ (P= Trin) x (P= Tox) | = P= Toy)? 


式 中 :7 为 日 间 平均 近 地 表 气 温 ;7 Pax A 了 分 
别 为 植被 光合 作用 最 低 、 最 高 和 最 适 气 温 , 也 可 以 利 
用 土地 利用 覆 被 数据 通过 查找 表 获 取 。 参 考 前 人 人 研 
jp 7 77 ,黑河 流域 植被 光合 作用 最 低 、 最 高 和 最 适 
气温 如 表 1 所 示 。 


表 1 各 生物 群 系 的 VPM 模型 参数 查找 表 
Tab.1 Reference table of parameters in VPM model 


for different biomes 


ES Emax T min Tias Top 
植被 类 型 /gC + (mol + APAR)-! PG rE / " 
常 绿 针 叶 林 0.42 -1 40 20 
落叶 阔 叶 林 0.42 -1 40 20 
FEHEM 0.42 1 48 31 
C3 草地 0.42 0 48 27 
湿地 0.42 si 40 20 

C3 作物 0.42 -1 48 30 
C4 作物 0. 63 -1 48 30 
人 造 建 筑 0.42 0 48 27 


基于 查找 表 与 土地 利用 覆 被 数据 ,生成 黑河 流 
域 最 大 光 能 利用 率 及 光合 作用 最 低 、 最 高 和 最 适 气 
温 的 空间 分 布 数据 ,使 用 聚合 法 实现 空间 分 辩 率 由 
30 m 降低 至 480 m, 

在 VPM 模型 中 ,水 分 胁迫 因子 W eaa p LSWI fi 
算 WARG) Hr TU : 

— 1 €LSWI 
sal plpewp. 

式 中 :由 下。 为 植物 生长 季 单 个 像 元 的 最 大 陆地 水 
分 指数 LSWI。 
3.2 ”模型 评估 

为 了 评 佑 本文 所 使 用 的 VPM 模型 的 性 能 与 总 
初级 生产 力 GPP 估算 的 精度 ,本 文采 用 了 2 种 方 
法 :直接 与 地 面 测量 CPP 相 比 ,与 MODIS CPP 产品 
(MOD17A2H) 进 行 比较 。 判 定 系数 (R ) 和 总 均 方 
根 误差 (RMSE ) 在 本 文中 被 用 于 量化 模型 估算 结 
gU, 


4 总 初级 生产 力 估算 结果 与 时 空格 局 


4.1 模型 估算 结果 评估 
通常 情况 下 , 通 量 监测 塔 的 测量 范围 有 限 ,一 般 
为 监测 塔 高 度 的 100 倍 ” 。 基 于 涡 动 相关 仪 高 度 ， 


(8) 
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考虑 像 元 和 站 点 测量 范围 的 共 配 准 误差 ,提取 VPM 
模型 与 MODIS GPP 产品 估算 的 站 点 GPP 值 。VPM 
模型 估算 和 地 面 实测 GPP 之 间 的 散 点 图 及 MODIS 
GPP 产品 和 地 面 实测 GPP 之 间 的 散 点 图 如 图 3 所 
示 。 

如 图 3a 所 示 ,VPM 模型 估算 的 GPP 和 地 面 实 
i| GPP 之 间 存 在 较 高 的 相关 性 ,判定 系数 RR 为 
0.849 95 ,总 均 方 根 误差 RMSE 为 1. 899 4 gC - 
m+ d7! , 拟 合 公式 为 y =0. 807 11x -0. 174 32, 
但 两 者 的 拟 合 曲 线 依旧 轻微 偏离 1: 1 参考 线 ,这 可 
能 是 地 表 植 被 的 空间 异 质 性 和 像 元 与 站 点 的 共 配 准 
误差 导致 的 。 单 个 像 元 30 m x30 m 范围 内 的 土地 
利用 覆 被 类 别 被 默认 为 是 单一 均匀 的 ,但 地 表 实 际 
情况 是 不 均匀 性 且 复 杂 的 。 同 时 , 通 量 贡 献 区 范围 
并 不 是 稳定 不 变 的 , 且 像 元 和 站 点 测量 范围 的 共 配 
准 误差 是 无 法 完全 消除 的 ,这 些 因 素 也 会 对 模型 估 
算 的 精度 造成 影响 。 

如 图 3b 所 示 ,MODIS GPP 产品 估算 的 CPP 和 
地 面 实测 CPP 之 间 相 关 性 较 弱 , 判定 系数 RR 为 
0.583 96, 总 均 方 根 误差 RMSE 为 4.418 3 eC - 
m^-*d' WEARY y 20.252 5x +0. 292 86, 与 
VPM 模型 估算 的 CPP 不 同 , MODIS GPP 产品 估算 
的 GPP 远 低 于 地 面 实测 CPP , 拟 合 曲线 显著 偏离 
1:1 参考 线 线 , 其 原因 可 能 有 2 种 :其 一 , MODIS 
GPP 产品 所 使 用 的 土地 利用 覆 被 类 别 的 最 大 光 能 
利用 率 su. 相对 较 低 ;其 二 ,MODIS CPP 算法 所 使 
用 的 土地 利用 覆 被 数据 MCDI2QI 在 复杂 地 表 条 件 
下 的 总 体 精 度 不 高 ,存在 错 分 的 情况 。 

因此 ,本文 使 用 VPM 模型 估算 的 总 初级 生产 力 
精度 要 高 于 MODIS CPP 产品 ,所 用 的 VPM 模型 要 
优 于 MODIS CPP 产品 所 用 的 MODIS PSN 模型 。 
4.2 黑河 流域 总 初级 生产 力 空间 分 布 

黑河 流域 生长 季 累 积 CPP 存在 着 明显 的 空间 
分 布 格局 。 沿 着 上 游 ,中 游 和 下 游 ,从 南 到 北 ,CPP 
的 分 布 具有 明显 的 梯度 ,如 图 4 所 示 。 

GPP 在 上 游 的 山区 南部 ,中 游 的 绿洲 区 域 和 下 
游 的 水 体 附近 的 河岸 区 域 相对 较 高 。 其 中 ,中游 以 
农作物 为 主要 植被 覆盖 的 绿洲 区 域 CPP 最 高 , 普 i 
KF 800 eC m”? o d ,上游 以 草地 为 主要 植被 覆 
盖 的 山地 区 域 GPP 次 之 ,大 多 小 于 600 gC +m”? - 
d ,下游 以 沙漠 和 戈壁 为 主要 土地 利用 类 型 的 地 区 
最 低 , 除 河岸 周边 外 大 面积 区 域 均 为 0 gC sm? - 
d ,这 种 空间 分 布 格局 主要 受 黑 河流 域 的 温度 和 
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图 3 VPM 模型 与 MODIS CPP 产品 的 估算 结果 评估 
Fig.3 Evaluation of estimated results from VPM model and MODIS CPP products 
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图 4 2015 年 黑河 流域 生长 季 累 积 CPP 空间 分 布 
Fig.4 Spatial distribution of accumulated CPP in Heihe 


River Basin during the growth season of 2015 


水 分 条 件 驱 动 : 上 游 山 地 地 区 ,年 平均 温度 为 -3.1 
~3.6 C ,降水 量 为 50 ~ 450 mm" ,植被 的 光合 作 
用 主要 受 低温 条 件 的 抑制 ;中 游 地 区 ,年 平均 气温 为 
7 -8.2 % ,人 工 灌溉 为 农作物 生长 提供 了 良好 的 水 
分 条 件 " ,环境 条 件 优 越 , 有 助 于 农作物 的 光合 作 
用 ;下 游 地 区 ,年 平均 降水 量 仅 为 40 ~ 60 mm ,但 年 
平均 潜在 莹 发 率 却 达到 2 300 ~3 700 mm" ,植被 
的 光合 作用 因 缺 水 而 受到 抑制 。 另 外 ,这 种 空间 分 
布 格局 还 受 士 地 利用 履 被 情况 的 影响 :由 于 中 游 绿 
洲 区 域 种 植 有 大 量 C4 农作物 (玉米 ) ,导致 该 区 域 


最 大 光 能 利用 率 相对 较 高 ,因此 总 初级 生产 力 远 高 
于 上 游 与 下 游 地 区 。 
4.3 黑河 流域 总 初级 生产 力 时间 动 态 

首先 ,研究 黑河 流域 全 境 植被 覆盖 区 域 (不 包 
括 裸 地 、 积 雪 与 冰 ) 的 平均 日 CPP 的 时 间 动 态 规律 。 
以 植被 叶 面 积 指数 表征 植被 生长 特征 ,如 图 5 所 示 。 
黑河 流域 全 境 植被 覆盖 区 域 的 平均 日 CPP 时 间 变 
化 趋势 与 黑河 流域 植被 生长 特征 基本 一 致 。 平 均 日 
GPP 最 初 随 着 时 间 的 推移 而 增加 ,并 在 7 月 20 H ~ 
7 月 27 H (201 ~209 qd) 达到 最 大 值 后 随 着 时 间 的 
推移 而 减少 。 该 变化 趋势 主要 是 受气 温水 分 .太阳 
辐射 和 植被 长 势 的 影响 而 导致 的 。 黑 河流 域 的 气 
温 、 降 水 量 与 河流 径流 量 在 7 月 与 8 月 间 达 到 一 年 
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. d! 
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图 5 植被 覆盖 区 域 的 平均 日 CPP 5 8 d & 
LAI 时 间 动 态 
Fig.5 Temporal dynamics of averaged daily GPP and 


8-day synthetic LAI in vegetation covered areas 
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图 6 黑河 流域 8 站 点 的 模拟 和 实测 GPP 时 间 动 态 


Fig.6 Temporal dynamics of simulated and measured GPP at 8 sites in Heihe River Basin 


中 的 最 高 值 '*” ,植被 叶 面积 在 7 月 间 达 到 娲 峰值， 
与 此 同时 ,太阳 辐射 在 6 月 达到 最 大 值 后 开始 逐 
渐 降 低 ,这些 因素 的 综合 影响 最 终 导致 了 黑河 流 
域 全 境 植被 覆盖 区 域 的 平均 日 CPP 的 倒 U 型 变化 
规律 。 

其 次 ,揭示 观测 站 站 点 CPP 的 时 间 动 态 规律 ， 
如 图 6 所 示 。 除 花 寨子 站 与 荒漠 站 外 ,所 有 观测 
站 的 日 CPP WER U 性 变化 规律 ,与 黑河 流域 全 
境 植被 覆盖 区 域 的 变化 规律 与 原因 基本 相同 。 而 
对 于 花 寨子 站 与 荒漠 站 ,由 于 其 下 热 面 为 荡 漠 与 
沙漠 ,植被 覆盖 率 极 低 ,因此 生长 季 的 CPP 总 体 趋 
势 维持 在 1 gC m-*. d-' 上 下 波动 ,季节 变化 不 
大 。 然 而 ,在 VPM 模型 中 ,对 于 低 植被 覆盖 的 裸 
地 ,由 于 其 最 大 光 能 利用 率 通 常 取 值 为 0 gC - 
(mol * APAR) ^! ,因此 GPP 的 估算 值 也 往往 等 于 
0 gCom“? + d^' , BOE WIE RAR AE RR TE 
盖 地 表 区 域 的 CPP 时 间 动 态 变化 规律 。 目 前 , 针 
对 低 植被 覆盖 地 表 区 域 的 最 大 光 能 利用 率 的 研究 
依旧 存在 较 大 不 足 。WANG 2579 曾 提 出 了 利用 最 
大 光 能 利用 率 与 增强 植被 指数 的 经 验 公 式 计 算 最 
大 光 能 利用 率 的 方法 ,但 该 方法 估算 的 最 大 光 能 
利用 率 误差 较 大 , 且 普 适 性 差 。 因 此 ,今后 对 低 植 
被 覆盖 地 表 区 域 的 最 大 光 能 利用 率 展 开 进一步 的 
深入 研究 是 有 必要 且 必 须 的 。 


5 讨论 


黑河 流域 地 处 西北 农 牧 过 渡 区 ,气候 类 型 为 
温带 大 陆 性 气候 ,植被 生态 系统 复杂 多 变 , 常 绿 针 
叶 林 落叶 阔 叶 林 ,高 寒 草 向 和 温带 草地 等 地 表 履 
盖 错 落 分 布 。 研 究 表明 ,VPM 在 常 绿 针 叶 林 、 热 带 
常 绿 林 、 落 叶 闪 时 林农 田 、 高山 草 多 和 温带 草原 
等 生产 系统 中 均 具 有 良好 的 模拟 效果 (R = 
0.64 ~0.94) ,因此 VPM 模型 在 黑河 流域 具有 和 较 
好 的 适用 性 。 本 文 利 用 VPM 模型 模拟 的 CPP 与 
通 量 站 的 观测 数据 进行 比较 分 析 (R =0.85) ,其 
结果 也 表明 了 VPM 模型 十 分 适合 用 于 在 黑河 流域 
GPP 模拟 。 

使 用 VPM 模型 进行 全 球 与 区 域 尺度 的 CPP 模 
拟 ,前 人 多 使 用 MODIS 产品 中 的 土地 利用 数据 
与 气象 站 点 插值 而 成 的 气象 数据 “| ,然而 前 者 精度 
在 具有 复杂 地 表 景 观 的 区 域内 较 差 ,后 者 精度 受气 
象 站 点 的 分 布 密度 影响 较 大 。 有 研究 表明 ,土地 利 
用 数据 对 地 表情 况 的 反映 随 着 分 辩 率 的 提高 而 更 加 
详细 准确 六 ,通过 陆 面 融 合资 料 可 获取 高 质量 与 格 
点 化 的 气象 要 素数 据 '" 。 黑 河流 域 地 表 和 覆盖 类 型 
复杂 , 且 气 象 站 点 分 布 较 东 南 沿海 相对 稀 下 ,使 用 高 
分 辨 率 的 土地 利用 数据 和 陆 面 融合 资料 可 以 大 大 提 
高 GPP 的 估算 精度 。 
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6 结论 


本 文 以 模拟 黑河 流域 2015 年 植被 生长 季 的 
GPP 为 例 ,时 间 与 空间 分 辨 率 分 别 设 为 8 d 与 500 m, 
使 用 MODIS 影像 数据 .ASTER GDEM 数据 、30 m 分 
状 率 土地 利用 覆 被 数据 和 中 国 区 域 地 面 气象 要 素 驱 
动 数据 ,驱动 VPM 模型 进行 CPP 估算 ,并 以 此 方法 
为 基础 ,揭示 黑河 流域 CPP 的 空间 分 布 与 时 间 动 态 
规律 。 人 研究 结果 表明 : 

(1) 本 文 VPM 模型 估算 的 结果 优 于 MODIS 
GPP 产品 的 预测 精度 ,以 大 沙龙 站 大 满 站 FEE 
站 .荒漠 站 混合 林 站 ,胡杨 林 站 .四 道 桥 站 和 湿地 站 
的 实测 通 量 数据 验证 ,其 中 ,后 者 的 判定 系数 仅 为 
0.583 96 , 157r 4812215 4.4183 gC em? - d^, 
而 前 者 的 判断 系数 达 0. 849 95, 提高 了 45.596 ,总 
均 方 根 误差 仅 达 1. 899 4 gC - m^ + d ,降低 了 
57.0% 。 

(2) 生 长 季 内 ,黑河 流域 上 游 ,中 游 和 下 游 的 
GPP 存在 着 明显 的 空间 分 布 梯度 ,主要 归 因 于 黑河 
流域 的 土地 利用 覆 被 情况 和 温度 与 水 分 环境 因素 ， 
其 中 ,中 游 地 区 的 水 热 条 件 最 佳 , 光 能 利用 率 高 , 生 
长 季 GPP 最 高 ,普遍 大 于 800 gC «m^ «d^; EYF 
地 区 的 光合 作用 受 低温 抑制 ,生长 季 CPP 相对 较 
低 ,大 多 小 于 600 gC «m^ - d^; FÜEREDCER T 5 
影响 而 地 表 缺 乏 植 被 覆盖 ,生长 季 GPP 最 低 , 除 河 
岸 周边 外 大 面积 区 域 均 为 0 gC - m^ +d, 

(3) 生 长 季 内 黑河 流域 全 境 与 局 部 有 植被 覆 
盖 的 区 域 (不 包括 大 片 沙漠 戈壁 ) 的 日 CPP 呈现 
先 增加 后 减少 的 倒 U 型 变化 趋势 ,该 变化 趋势 主 
要 是 受 植 被 长 势 .气温 .水 分 和 太阳 辐射 的 影响 而 
导致 的 ,其 中 ,全 境 植被 覆盖 区 域 的 日 CPP 在 7 月 
下 旬 (7 月 20 日 ~7 月 27 日 ) 达 到 最 高 值 , 最 高 值 
约 为 5.65 gC- m””. d”, 与 植被 生长 特征 基本 保 
持 一 致 。 

(4) 生 长 季 内 黑河 流域 境内 的 沙漠 、 议 漠 \ 裸 地 
等 植被 覆盖 率 极 低 的 土地 日 CPP 值 在 基本 稳定 中 
上 下 波动 ,稳定 值 大 致 处 于 1 gC +m? + do" 附近。 
光 能 利用 率 模型 在 估算 该 类 土地 利用 类 型 时 ,常常 
默认 为 0 gC * (mol + APAR) ^ ,不 能 反应 其 时 间 动 
态 变 化 。 

本 文成 功 的 估算 了 2015 年 黑河 流域 的 总 初级 
生产 力 ,然而 研究 沿 有 不 足 之 处 存在 :(1) 研究 时 间 
段 较 短 ,不 能 全 面 的 揭示 研究 区 CPP 的 年 际 时 空 变 


化 ;(2) 空间 分 辩 率 较 低 ,空间 蜡 质 性 对 估算 精度 的 
影响 依旧 较 大 ; (3) 仅 对 C3/C4 植被 的 最 大 光 能 利 
用 率 进 行 区 分 ,未 能 针对 不 同 的 生物 群落 采用 不 同 
的 最 大 光 能 利用 率 。 因 此 今后 有 必要 进一步 扩大 研 
究 时 间 段 ,提高 空间 分 辩 率 ,人 研究 最 大 光 能 利用 率 与 
生物 群落 的 联系 ,从 而 对 研究 方法 作 进一步 改进 。 
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Estimating gross primary production in the Heihe River Basin 


from multiple data sources 


TONG Zhi-hui', XIONG Zhu-guo', SUN Rui^^, LIU Xiang-tong', WANG Dong-dong! 
(1 College of Surveying and Mapping , East China University of Technology , Nanchang 330013 , Jiangxi , China ; 
2 State Key Laboratory of Remote Sensing Science , Beijing Normal University , Beijing 100875 , China; 

3 Beijing Engineering Research Center of Global Land Remote Sensing Product , Beijing 100875 , China ) 


Abstract: The initial matter and energy entering terrestrial ecosystem are determined by the gross primary produc- 
tivity (GPP). Hovexer ,the GPP measurements acquired by the MODIS GPP products are not sufficient to accurate- 
ly reflect the distribution of ecosystem matter and energy in the Heihe River Basin in Qinghai Province ,Gansu Prov- 
ince and Inner Mongolia, China with complex surface coverage. Therefore , based on MODIS image data, ASTER 
GDEM data,land cover data with a spatial resolution of 30 m, and the China Meteorological Forcing Dataset , the 
VPM model was derived to simulate the gross primary productivity of the Heihe River Basin from May to October 
2015 with a spatial resolution of 500 m and a temporal resolution of 8 days. Based on simulations using the VPM 
model,the spatial and temporal patterns of GPP in the Heihe River Basin during the growing season were deter- 
mined. The results of the study indicate that the accuracy of the results estimated using the VPM model was higher 
than that of the MODIS GPP products. The judgment coefficient increased by 45. 596 , while the total root mean 
square error reduced by 57.096 . The study results also demonstrate that the GPP accumulation during the growing 
season in the Heihe River Basin exhibited a significant spatial distribution gradient pattern , which can be described 
as the highest in the middle reaches ,the next highest in the upper reaches , and the lowest in the lower reaches. In 
addition ,the daily GPP of the whole and partial vegetation-covered areas in the Heihe River Basin first increased 
and then decreased in an inverted U-shape. The daily GPP values of all vegetation-covered areas reached the maxi- 
mum in late July. The daily GPP values of those ground areas with very low vegetation coverage fluctuated up and 
down in their basic stability , while the stable value was approximately 1 gC * m! * d^'. 

Key words: VPM model; multi-source data; gross primary productivity; spatial and temporal patterns; Heihe 


River Basin 


